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^1^' We apply the fractional exclusion statistics (FES) to the quantum hadrodynam- 

j-^. ics model (QHM). The FES typically describes in an effective way the remanent 

\ interaction between the constituent particles. We determine the phenomenological 

' parameters of the QHM with FES from the conditions at the saturation point of the 

cn \ 

OO ■ symmetric nuclear matter, independent of the exclusion statistics parameter, a. We 

; 

, calculate the relevant physical quantities of the system, e.g. the effective mass of the 

00 ' nucleons and its dependence on the baryon density, temperature, and a. We also 

o : 

^ ' calculate the isothermal curves (pressure vs. baryonic density) for different temper- 

^ . atures and as. We observe that in the region of phase transition the isotherms do 

H : 

5t , iiot depend much on the exclusion statistics parameter, but a strong dependence on 

a is observed at higher densities or at higher temperatures. 
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INTRODUCTION 



The relativistic mean-field model (RMF) is a general framework for the relativistic nuclear 

ndnri 

many-body problem based on hadronic degrees of freedom [l|, 12, y, |J] . The RMF model has 
successfully been used to describe many nuclear phenomena [51]. In the models of nuclear 
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matter studied so far, the statistics of particles involved was not a subject of discussions. 
The baryons are fermions and mesons are bosons. However, in the hot and dense nuclear 
medium the statistics of quasi-baryons may change. Therefore, the main purpose of this 
paper is to investigate the influence of the statistics of particles on the thermodynamic 
properties of the isospin symmetric nuclear matter. For this reason, we apply the concept 
of the fractional exclusion statistics to the RMF of Refs. s], 4|. In FES, the statistics of the 
ideal particles is deformed by the additional parameter a, and so it is posible to fit more of 
the known properties of nuclear matter. 

The fractional exclusion statistics, introduced by Haldane in Ref. jsl has received very 

las been applied to many models of interacting 

m 



much attention since its discovery and 



systems (see 
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or example Refs. [71, | 
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lany models of int( 
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28|). Several authors have a. 



manifestation of FES 
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so discussed the microscopic reason for the 
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33|. Iguchi and Sutherland 



33(1 showed that liquids of particles in three dimensions, interacting through long-range 
forces exhibit the nature of quantum liquids with FES, the characteristics of the FES being 
determined by the interaction. Murthy and Shankar 20|] analysed a system of fermions 
in the Calogero-Sutherland model and showed that by redistributing in a certain way the 
interaction energy in the system among the single-particle states, one obtains a gas which 
is described by FES. 

Although the concept received so much attention and has been applied in general to 
many types of systems, in Ref. 3^ it was proven that it has basic inconsistencies. These 
inconsistencies have been corrected in the same paper, by introducing mutual exclusion 
statistics parameters which are proportional to the dimension of the subspace on which they 



act. In a related paper it was also shown that the FES is not an exceptional statistics, 
manifesting only in certain exotic systems, but it is generally present in interacting particle 
systems. Concretely, a non-ideal system of interacting particles may be described as an 
"ideal" FES system, in which the exclusion statistics parameters are determined by the 
interaction Hamiltonian. 

The structure of the article is the following. In Section [Til we briefly describe basic 
ingredients of the relativistic mean-field model in fractional exclusion statistics approach 
and the methodology to evaluate the model parameters. The thermodynamic results are 
discussed in Section III CI The main conclusions are summarized in the final Section IIIII 
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Throughout the paper we use the natural units system, fi = c = ks = I- 

II. FRACTIONAL EXCLUSION STATISTICS APPLIED TO QUANTUM 

HADRODYNAMICS 

A. The relativistic mean-field model 

As already mentioned, RMF is a relativistic quantum theory model for the nuclear many- 
body system. The motivation to introduce the RMF was the observation of large Lorentz 
scalar and four-vector components in the nucleon-nucleon interaction. 

To introduce the notations and write the basic equations, we follow the treatment from 
In the model we have the barionic field, ip, that describes the nucleons, p and n, the 
neutral scalar meson field, </>, for the meson a, and the neutral vector meson field, V"^. 
We assume that the nuclear scalar meson couples to the scalar density of barions through 
gs'4'ip(f) ^^'^ that the neutral vector meson couples to the conserved barion current through 

We also assume that the baryon density is high, so that we can use mean field theory 
(MET) and replace the meson field operators by their expectation values, which are the 
classical fields (0) = 0o and (V^) = S^qVq. For a static, uniform system, the quantities 
00 and Vo are constants, independent of x^. With these simplifications (see Refs j^, 4] for 
details), we have 0o = {'4^'^)9s/f^1 = Ps9s/f^1 and Vq = {ip'^ip)gi,/ml = psgv/fnl, with 
rris and rriy the masses of the scalar and vector mesons, respectively, and the mean field 
Lagrangian is written as 

£mft = ^ - g^n^'Vo - (M - gs<Po)] ^ - \rnl4>l + ]^mlVl (1) 

Applying the Lagrange equations to £mft one obtains the equation for the baryonic Dirac 
field, 

[il,d^ - ^7.7Vo - (M - gs<i>o)]ip = 0, (2) 
which, when solving for plane-waves, gives the single-particle spectrum 



= E*{k)+g,V^, (3) 
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and 

e{k) = E*{k)-g.,Vo (4) 

with M* being the effective mass, 

M*=M- gs<j)o. (5) 

The baryonic number, 

73^ rO _ / ^3^,/,t 



B = J^d'^B'' = J^d':>^'^^i/j = VpB, (6) 

is a constant of motion, so it is conserved. (Throughout the paper, we shall denote the 
volume of the hadronic system by ^-without any index.) 

From the energy-momentum tensor calculated with ([T]) one can derive the expressions for 
the pressure and energy density, 

p = ^^^{-la ■ V)^ + ^mlV^^ - ^m^cf^l, (7) 

£ = ^t [_,^ . V + + g^Vo] - ^m'X + ^"^'0o- (») 
The total Hamiltonian of the system is 

H = rf^xToo = d^^S. (9) 
In the second quantization, the baryon number and the effective Hamiltonian operators 

are 

B = ^(4,AkA-5Li?kA), (10) 
k,A 

i^MFT = E Vk2 + M*2(4,AkA + Bl,B^,) + g^VoB - V (^^m'X ' ^^s'/'o) , (H) 

kj A 

where Al^^^ and are the creation and anihilation operators for the baryon state of mo- 
mentum k and species A (A = n or p, of spin up or spin down), whereas B^^^ and i?kA are 
the corresponding antiparticle operators. 



B. Thermodynamics 



We shall consider in this paper the FES in its simplest form, which is the diagonal form. 
We introduce the spin-isospin degeneracy factor, 7, which takes the value 2 for neutron 
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matter (two spin projections and one isospin projection) and 4 for nuclear matter (two spin 
and two isospin projections). The net baryon number i? is a conserved quantity, so the 
chemical potential, /i, is associated to it. In the grand canonical ensemble the partition 
function is written as 

Z = Tr (^e-^^^^) . (12) 

Plugging Eqs. (fTOj) and (fTT]) into (fT2!) and maximizing the partition function, we obtain in 
the standard way the thermodynamic quantities, 

PB = JjMT) - nkiT)]d% (13) 

^ - 2|^^-i^^-^*)^+3(2^L^^k^M^^ (1^) 
where nk{T) and nk(T) are the baryon and antibaryon mean occupation numbers of the 
single particle level k. We assume that all the baryons have the same mass, and therefore 
their spectra and single-particle level populations are identical. 

If the baryons are ideal fermions, then nfc(T) and nk{T) are simply the Fermi functions, 

where i/ = /i - g^Vo = /i - glpB/ml. 



In general, in a 
but by the formula 



ES system the mean occupation numbers are not given by Eq. f|T6l) . 



29 



30 



34| 

nk,a = [w{i) + a\-\ (17) 



where w is a solution of the equation 

w"(0[l + u;(0]'"" = ^ = e(^-^)/^ (18) 

(note that we added the subscript a to specify the statistics). In this case the formalism 
remains the same, except that we work with the more general levels population (1171) . instead 

of (USD. 

The equilibrium values of the fields 0o cind Vq, and therefore of M*, are determined by 
maximizing the grandpartition function, Z, with respect to 0, which gives js, ^ 

m2(27r)3yfi3 Vk2 + M*2 ^ ' 
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that has to be solved for M* . 

The phenomenological constants of the model, expressed as = g'^M'^/m'j. and = 
g^M"^ /ml are determined from the conditions at the saturation point: at T = and pb = Po, 
the binding energy, Sb = [£{pb) — PbM]/pb attains its minimum (i.e. equilibrium) value, 
Sbo- For the nuclear matter we take [4] 

£bo ^ -16 MeV, Po = 0.16 fm'l (20) 

Therefore we rewrite Eqs. ( fT3l) - (|T5l) . and ( |T9|) at T = 0: 

7 f .31. 7 ,3 



[ZTTj-^a Jk<kp bn'^a 

S = /^p| + ^(M - M*f + -P— [ d'kE*ik), (22) 
2mfj 2gj (27r)-^a Jk<kp 

2 ^2 ^. , 1,2 



9v 2 '^s/i^j T\/r*\2 , 7 /■ ,3i_ 



= M-M*-^-4^/ d^k — — , (24) 
m2 (27r)3a A<fe^ E*(A;)' ^ ^ 

where /ci? is the Fermi wavevector. If we minimize Sb for the usual Fermi statistics, a = 1, 

and with the expression for S given by fl22|) . we obtain ~ 330 and ~ 249. If we keep 

Cs and Ct, fixed at these values and plot Sf, for the nuclear matter, as a function of p^, for 

different as (see Fig. [T]a), we observe that the binding energy decreases and its minimum 

moves towards higher baryonic densities; so the characteristics of the nuclear matter change 

significantly with the exclusion statistics. This situation have been analysed in Ref. 35 1. 

Since the exclusion statistics is a microscopic characteristic of matter, which is employed 

(typically) to describe the eventual remanent interaction in the MFT in an effective, statis- 



tical manner [28|, in this paper we assume fixed the experimentally observable quantities, 
Sbo and Pq, from which we determine the microscopic parameters Cg and as functions of 
a. 

Therefore we determine the parameters Cg and C^, so that the minimum of the function 
^b{pb), at T = 0, does not change with a. The condition of minimum for Sb reads 

d (S-pbM\ 1 dS 



PO ^0 



^ = 0. (25) 



dpB \ Pb ) Po dpB 
and calculating the total derivative of £ with respect to pb from Eq. ( 122|) . we obtain an 
equation for 

C?, Sm + M -E*(kpn) 
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where kpo = [(27r)^a;po/7]^''^- Using Eq. ( 126|) we eliminate Cy from the expression of £ ( 122 

£ko + M-E*{kFo) , M\M-M*f 
^ — Pb H 



+ 



7 



2 ' 2a 



M* 



(27) 



where fci;' is given by Eq. ( 12T1) as a function of and we evaluated analytically the integral 
over k. From the effective mass equation (l24l) . we get Cg as a function of M*, 

-1 



47r2a M^{M-M* 
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kFJkl + M*^ - M*^ In 



(28) 



which we insert into (1271) to obtain a self-consistent equation for M*: 

kFQE*{kFQ) - M*'^\n 



£o ^ PoE*{kFo) ^M*{M-M* 
2 2 



,2, kFO + E*{kFo) 



M* 



+ 



7 



r*4 1 fi-Fo 



kFoE*{kFo)i2k'pQ + M*') - M*Mn 



E*{k 



FO) 



M* 



(29) 



With M* calculated from the equation above, we go backwards and calculate Cs(a) and 
C^{a) from Eqs. (EHD and (EE]). 



C. Results 

In Fig. [1] (b) we show the binding energy of the particles in the nuclear matter for 
a = 1, 0.5, and 0.1. We see there that all the curves £h,a have their minima in the same 
point and that the binding energy start to depend on a roughly at Pb > Po- 

In Fig. [2] we show the dependence of Cg and on a, whereas in the inset we plot the 
ratio between the effective and the free-particle masses, M* /M. As a increases to 1, Cg{a) 
and C^(«) decrease monotonically to their values corresponding to the Fermi statistics, 
C2(l) ^ 330 and ^^(l) ^ 249. 

To calculate Cs and at a = we have to calculate first the asymptotic expression of 
M*, for a ^ 1. We notice from Eq. fl29l) that if M*^q is finite, then the second and the 
third terms at the r.h.s. of Eq. ( 1291) diverge. But since the l.h.s. is a constant and the first 
term at the r.h.s. converges to —PqM*^q/2, we conclude that Eq. ( l29i) cannot hold in such 
therefore M* must converge to zero as a — > 0. 

We can also notice by direct inspection of Eq. (!29l) that at small values of a, M* must 
have an asymptotic behavior M* ~ a^/'^, which is similar to the a dependence oi kFo- If we 
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denote M* = CM<y^^^ and kpo = cpa^^^ = [(27r)^apo/7]^^^ and plug them into Eq. ( 129|) . we 
obtain an equation for cm- 



I - ^ • {c.y4T4;-c-.ln '-"f; + ^h . (30) 
Using (!30i) we obtain 



^— and (0) = — 
Pb ^0 



C,^(0) = ^ and Ct{0) = —. (31) 



Note that when a converges to zero, M*/M converges also to zero, whereas Cg and C„ 
converge to finite values. 

In Fig. [3] we plot the ratio M*/M for nuclear matter (a) and neutron matter (b) as a 
function of temperature. In both plots the three upper curves correspond to zero baryon 
density {pB = 0), whereas the lower curves correspond to the barion density po- We notice 
that even though the values of M*(T) are different for different as, the general behavior is 
similar, i.e. for pb = 0, M* is a monotonically decreasing function of T, whereas for finite 
baryon densities, M* first increases and then decreases with T. 

In Fig. m we plot the isotherms p{pb), for three values of T (T = 10, 15, and 100 
MeV) and a = 1, 0.5, and 0.1. We notice that at low temperatures and baryon densities 
the isotherms depend very little on the exclusion statistics parameter, whereas at high 
temperatures the isotherms become very sensitive to a. This is due to the fact that the 
phenomenological constants Cg and are recalculated for each a, so that the saturation 
point is independent of the exclusion statistics. On the contrary, if the Hamiltonian of the 
system would be independent of a, i.e. we would keep Cs and the same for all as, then 
the high temperature results would coincide (statistics has smaller and smaller infiuence on 
the thermodynamic results as the temperature increases), whereas at low temperatures the 



isotherms would have a strong dependence on a 35|. 



III. CONCLUSIONS 



We applied the fractional exclusion statistics 61 to the relativistic mean field model of 

n ri 

Walecka and Serot [31, Fixing the observable (i.e. physical) parameters of the system, 
which are the binding energy, S^, and the baryonic density, Pq, we determine the model's 
phenomenological mean field parameters, which are Cs and Cy (see Fig. [1]). Therefore in 
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our model, Cg and are functions of the exclusion statistics parameter, a, as depicted in 
Fig. [21 Using Cs{a) and Cy{a) we calculated the effective masses of the nucleons in nuclear 
matter (Fig. [3] a) and neutron matter (Fig. [3]b). 

Using further the effective mass in the expression ( fT5i) for pressure, we draw in Fig. H] 
the isotherms for nuclear and neutron matters. We observed that at low temperatures and 
baryonic densities the isotherms depend very weakly on a, whereas at high temperatures the 
dependence on a is strong (Fig. H]). This dependence is different from the model with fixed 
Cs and Cy, in which case the a dependence of the isotherms is strong at low temperatures 
and weaker and weaker as the temperature increases. 

The phase transition appears in systems of any FES, but depends on the exclusion statis- 
tics parameter. 
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FIG. 1: The energy of the nuclear matter, at zero temperature: (a) the constants Cs and Cv are 
fixed at Cs = 330 and = 249, which corresponds to the binding energy E'to = — 16 MeV and 
the baryon density po = 0.16 fm~^ for Fermi statistics (a = 1); (b) the constants Cg and C„ are 
calculated for each value of a, so that the observables E'feo ^iid po remain the same. 
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FIG. 2: The constants Cg and C^, as functions of the exclusion statistics parameter, a, when the 
binding energy E^o = — 16 MeV and baryon density po = 0.16 fm"'^ are independent of a. In the 
inset is plotted the effective mass vs. a, under the same conditions. 
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FIG. 3: Effective mass, M*, as a function of the temperature, for a = 1, 0.5, 0.1 and for two values 
of pb'- pb = (the upper three curves in each graph) and pb = Po = 0.16 fm~'^ (the lower three 
curves in each graph). The constants Cs and are determined so that the binding energy and 
the net baryonic density at T = are EhQ = —16 MeV and Pb = Po = 0-16 fm~'^ for each a. The 
curves in (a) and (b) correspond to nuclear (7 = 4) and neutron (7 = 2) matter respectively. 




FIG. 4: Isotherms, p{pb), for three temperatures, T=10, 15, 100 MeV (the isotherms corresponding 
to T = 100 MeV are given in the insets), and a = 1 (solid line), 0.5 (dash line), 0.1 (dot line). Up: 
nuclear matter, 7 = 4; low: neutron matter, 7 = 2. 



